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Abstract 

We compute one of the second order transport coefficients arising from the chi¬ 
ral anomaly in a high temperature weakly coupled regime of quark-gluon plasma. 
This transport coefficient is responsible for the CP-odd current that is proportional 
to the time derivative of the magnetic field, and can be considered as a first cor¬ 
rection to the chiral magnetic conductivity at finite, small frequency. We observe 
that this transport coefficient has a non-analytic dependence on the coupling as 
~ 1/{g‘^\og{l/g)) at weak coupling regime, which necessitates a re-summation of 
infinite ladder diagrams with leading pinch singularities to get a correct leading log 
result: a feature quite similar to that one finds in the computation of electric con¬ 
ductivity. We formulate and solve the relevant CP-odd Schwinger-Dyson equation 
in real-time perturbation theory that reduces to a coupled set of second order dif¬ 
ferential equations at leading log order. Our result for this second order transport 
coefficient indicates that chiral magnetic current has some resistance to the time 
change of magnetic field, which may be called “chiral induction effect”. We also 
discuss the case of color current induced by color magnetic field. 


*e-mail: Emiadeo.j@gmail.com 
le-mail: hyee@uic.edu 



1 Introduction 


The chiral anomaly is an intriguing quantum mechanical phenomenon arising from an 
interplay between charge and chirality of massless particles such as chiral fermions. It 
has recently been appreciated that the chiral anomaly may induce interesting parity-odd 
transport phenomena in the plasmas of such particles [1, 2, 3, 4, 5]: at the lowest order in 
derivative expansion of hydrodynamics (that is, at the first order) it has been shown that 
the second law of thermodynamics dictates the existence of such phenomena [6]: the chiral 
magnetic effect [3] and the chiral vortical effect [7, 8]. Moreover, the magnitudes of these 
transport phenomena in the static, homogeneous limit are hxed by underlying anomaly 
coefficients and are not renormalized by interactions. This universality has been conhrmed 
explicitly in both weak [9, 10, 11, 12] and strong coupling [13, 14, 15, 16, 17] computations, 
and there are also evidences in favor of them in lattice simulations [18, 19, 20, 21, 22]. 
Recent results from heavy-ion experiments at RHIC [23, 24, 25, 26, 27] and LHC [28] 
seem consistent with the predictions from chiral magnetic and vortical effects [29, 30] 
(as well as chiral magnetic wave [31, 32, 33, 34, 35]), and quite interestingly, there is 
a successful experimental test of chiral magnetic effect in Dirac/Weyl semimetals which 
feature chiral fermionic excitations [36] (the spin degree of freedom in this case arises 
from an internal degeneracy, in other words, it is a pseudo-spin). Therefore, the existence 
and the magnitudes of these transport phenomena at lowest order (i.e. hrst order) in 
derivative expansion seem by now quite robust. One can also generalize them to all even 
higher space-time dimensions than four [37, 38, 39, 40]. 

As we go beyond the lowest order in derivatives, the possible anomaly induced trans¬ 
port phenomena become numerous: in four dimensions there are thirteen possible second 
order anomalous transport coefficients in the current and energy-momentum tensor in a 
conformal plasma [37] (and more in non-conformal plasma [41]), while the second law of 
thermodynamics seems to constrain only eight combinations of them [37]. Some of these 
constraints have been confirmed in a strong coupling computation [42]. The interesting 
fact is that the values of these anomalous second order transport coefficients, although 
they are proportional to anomaly coefficients, do depend on the dynamics of the micro¬ 
scopic theory up to the mentioned constraints, so computing them in weak and strong 
coupling regimes is a non-trivial, but worthwhile task in any theoretical model. 

The purpose of this work is to take a small step in computing these second order 
anomalous transport coefficients in weakly coupled gauge theories, having in mind QCD 
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and electro-weak theory. Our current study will be based on diagrammatic techniques, and 
we hope to address a similar computation in chiral kinetic theory framework [43, 44, 45, 46, 
47] in a separate work. We will show that one particular second order anomalous transport 
coefficient in the charge current has a non-analytic dependence on the coupling constant, 
~ 1 /\og{l/g), which is similar to that one hnds in the shear viscosity and electric 
conductivity [48, 49] (and also in the chiral electric separation conductivity [50])*. This 
transport coefficient appears in the second derivative correction to the current constitutive 
relation as, 

~ ^5 , ( 1 . 1 ) 

where is the electric field strength in a local fluid rest frame dehned by 

and we followed the notation introduced in Ref. [37] to denote the transport coefficient 
,^ 5 . Using the Bianchi identity, one can replace e^''°‘^Uy'DaEp with with = 

being the magnetic field strength in the local rest frame, which means that .^5 
can be viewed as a hrst correction to the static chiral magnetic effect at hnite frequency. 
More explicitly, it appears in the anomalous part of the current density as 

^ ^ dB 

J = ■ ■ ■ + a^B + ^5-^ + ■ • • ; (1-2) 

where is the topologically protected value of chiral magnetic conductivity at zero 
frequency. .^5 is parity (P) and charge conjugation-parity (CP) odd, so must arise from 
chiral anomaly. 

As was observed hrst by Jeon [51], in diagrammatic language, the non-analytic behav¬ 
ior in the coupling dependence is signaled by the presence of pinch singularities in multi¬ 
loop ladder diagrams of two point correlation functions, which necessitates a resummation 
of all ladder graphs by solving a Schwinger-Dyson type equation to get a leading log result. 
Previously, it was observed in Ref. [12] that the zero frequency-momentum limit of the P- 
odd part of the 1 -loop correlation function does not have pinch singularity, reproducing 
the correct static value of chiral magnetic conductivity. We hrst motivate our study by 
observing an appearance of pinch singularity in the P-odd part of 1-loop diagram at hrst 
order in frequency, which enters in a Kubo formula for ,^ 5 . Following intuitions from the 
computation of electric conductivity [52, 53, 54], we then identify multi-loop ladder graphs 

*Our definition of transport coefficients does not include a trivial factor from the definition of 
electromagnetic current which is e times of the fermion number current. Therefore, all quantities in 
our work are defined with the fermion number current. For example, the electric conductivity will be 
~ l/(e^ log(l/e)) and the chiral magnetic conductivity at zero frequency for a single right-handed Weyl 
fermion is 
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whose P-odd parts contain a chain of pinch singularities that have to be re-summed to 
get a correct leading log result for ,^ 5 . The emerging Schwinger-Dyson equation is more 
difficult to solve than that in the electric conductivity, because we need to keep hnite 
external momentum k (up to hrst order in k) to extract P-odd part of the correlation 
function. In section 4, we prove an important fact that all k dependence in the denomi¬ 
nators of pinching propagators do not contribute to the P-odd part of our interest up to 
first order in k, allowing us to neglect them in the denominators of pinching propagators. 
The necessary k dependence for P-odd correlation functions arises only from the spinor 
projection part of the fermion propagators. With this important simplihcation, we are 
able to reduce the leading log part of the P-odd Schwinger-Dyson equation into a coupled 
set of second order differential equations, which can be solved numerically. Along the 
way, we develop and use the sum rule for the P-odd part of Hard Thermal Loop (HTL) 
photon spectral density, which is summarized in the Appendix 1 . 

For most part of our presentation, we will consider a single species of Weyl fermion 
in quantum electrodynamics (QED) for simplicity, and a generalization to hnite number 
of species of Weyl and Dirac fermions as well as to a non-abelian SU (A^c) gauge theory 
is trivial at our leading log order. We will describe this generalization in our discussion 
section at the end. Our results are summarized as follows: for QED with a single right- 
handed Weyl fermion, we have 

^ —3.006 /i 

^ e4 log(l/e)T ■ 


For 2-havor massless QCD (Aj, = 3) with = 3/2 and Qd = —1/3, our result is 


tQCD 

S5 


(Qn + Ql) 


-3.6 /iA 
^Dog(l/^) T 


-2.003 /iA 
^Dog(l/^) T ■ 


(1.4) 


where /ta is an axial chemical potential. The sign of .^5 compared to the zero-frequency 
value o'^(= ^ for QED) is a meaningful dynamical result. A relative negative sign 
between the two means that the chiral magnetic current has some resistance to the change 
of the magnetic held. We may call this “chiral induction ehect”. 


2 Pinch singularity in P-odd part at 1-loop 

In this section, let us motivate our work by observing an appearance of pinch singularity 
in the P-odd part of 1-loop diagram at hrst order in frequency u. It will also serve to 
hx our notations and conventions. For simplicity, we will consider the case of single Weyl 
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fermion species of unit charge in QED plasma at finite equilibrium temperature T, as 
the generalization to multi flavors or non-abelian gauge groups is simple (which will be 
summarized towards the end of the paper). Throughout our analysis, we will use the real¬ 
time Schwinger-Keldysh formalism in “ra”-basis to compute the retarded current-current 
correlation function that contains the chiral magnetic conductivity (T^{k) in its P-odd 
part, 

= ia^{k)e^^^k ^, (2.5) 

where italic letters run over the three spatial dimensions, and k = (A;°, k) = k) is an 
external four momentum. Note that in ra-basis, the retarded two-point function is equal 
to 

(r{k).P{-k))^ = {-i){J‘{k)Ji{-k))sK = {-i)GlJk) , (2,6) 

where the subscript SK in the second term emphasizes that it is computed in the Schwinger- 
Keldysh path integral with = l/2(Ji -|- J 2 ) and = Ji — J 2 (1 and 2 denote the two 
time contours in the Schwinger-Keldysh formalism). We follow the notations in Ref.[12] 
for consistencyb Explicitly, 

Jjf = , (2.7) 

with (j^ = ( 12 x 2 , <7) in terms of two component Weyl spinor held Therefore, the task 
is to compute the P-odd part (or anti-symmetric part in i, J indices) of the (ra)-correlator 
^(ra)(^) small frequency-momentum. 

The zero frequency-momentum limit of <J-^{k) has been shown to be universally 

lim lim aJk) = , ( 2 . 8 ) 

and in particular, there appears no pinch singularity in this limit as shown in Ref. [12]. 
Our ,^5 appears in hrst order expansion inuj = kP (while still taking zero momentum limit 
k ^ 0), 

\ima^{k) = ^ - . (2.9) 

However, since the P-odd part of contains a linear term in k in dehning (T^{k), 

= -’’xikV’V . ( 2 . 10 ) 

we have to keep k dependence in up to hrst order in k, and then take fc —)■ 0 

limit of the P-odd coefficient o'^(/c). This essentially means that we need to keep hnite 

Tn literature it is often chosen to denote the retarded function by G'(ro), which we hnd confusing. 
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k in the Schwinger-Dyson equation for ladder resummation, which is in contrast to the 
case of electric conductivity where one can put fc = 0 from the very outset which greatly 
simplihes the analysis. Despite this difficulty, we will be able to solve the Schwinger-Dyson 
equation for the P-odd part of and extract the coefficient .^ 5 . At 1-loop there 



k 

<- 


Figure 1; Diagrams responsible for the retarded response at one loop in the ra-basis. 


are two Feynman diagrams for in real-time formalism as depicted in Fig.l. The 

fermion propagators are given by 

i 


Sra{p) = 


- s|p| +iC/2 


S = ± 


Sarip) = 


s=± 


s | p | -</2 


pu 


Vs{p) , 


Vs{p), 


Srr{p) = (^2 “ j P(P) ’ 

where n±(p°) = 1 ), and the spectral density p is 

c 


d(p) = 


t^{p^-t\p\f + {Cm- 


Vtip), 


( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 


and we introduce the damping rate C ~ log(l/p)T in the propagators, which will be 
needed to regularize possible pinch singularities, that is essential to have a non-analytic 
dependence on the coupling constant. It is important to observe a thermal relation 


SrAP) = ( § “ n + (p°) ) (Sra(p) “ S„r(p)) , 


( 2 . 15 ) 


which plays a central role in our analysis. The projection operators V± {p) are dehned as. 


'P±{p) = - 1 ± 


cr ■ p 

\P\ 


= T 


(y ■p± 

2|P| ’ 


(2.16) 
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where = ( 1 , —cr) and = (±|p|,p). Our metric convention is (—,+,+,+). The 
operators V±{p) project onto particle and anti-particle states respectively with given 
momentum p, and the {s,t) summation in the above physically represents distinctive 
contributions from particles and anti-particles. 

The 1-loop expression for from the two Feynman diagrams is 

(“1) j [<^"Sra{P + k)a^Srr{p) + (j"Srr{p + k)a^ Sar{p)\ ■ (2.17) 

where (—1) in front comes from fermion statistics. In Ref. [12] it was shown that after 
extracting e^^^k^ for the P-odd part, the limit fc —)■ 0 and ca —)■ 0 commutes and produces 
the correct result /i/ 47 r^, without featuring pinch singularity. However, we will see that 
ca —)■ 0 limit hides the pinch singularity appearing at hrst order in u for ,^ 5 . Using the 
thermal relation (2.15) to replace Srr with {Sra — Sar), we have several combinations of 
Sra and Sar- From the well-known fact (see Ref. [51]) that the pinch singularity appears 
only from the pair of Sra and Sar sharing a same momentum^, let us select only terms 
that potentially contain pinch singularity, which results in 

= y {n+{p^+ u)-n+{p^)) tr [a"Sra{p + k)a^Sar{p)] ■ (2.18) 

It is clear from this expression that a; —)■ 0 limit does not produce a pinch singularity, 
because (n+(p° -|- cj) — n+(p°)) ~ {dn+{p^)/dp^)u + 0(u‘^) already gives a linear factor in 
u. Moreover, one can put a; —)■ 0 in the rest expression as we are only interested in the 
linear term in u for ,^ 5 . In computing the above using (2.11) and (2.12) for Sra and Sar, 
let us recall that the chosen s = ± representing particle or anti-particle from Sra must be 
the same s chosen in Sar to have a pinch singularity in their product. Therefore, we have 

G;';“(A;) ^ + k)a^Sar{p)] 

^ f d^P f dn+{p^)'\ f tT[cT^Vs{p + k)a^Vs{p)] 

^ j \ dp^ ) ^^{p^ — s\p + k\ + iC^/2){p^ — s\p\ — 1 (^/ 2 ) ' 

(2.19) 


To identify the P-odd structure containing e^^^'k!' in small fc —)■ 0 limit, we hrst note 
that there are two possible sources of fc-dependence: one is from the denominator and the 

Mhis is because Sra {Sar) bas particle poles slightly below (above) the real axis by an amount ±i^/2, 
so that the residue of their product contains 1/^ factor which is the (regularized) pinch singularity. 
This also means that s = ± in (2.11) and (2.12) must be common in Sra and Sar pair causing a pinch 
singularity. 
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other is from the projection operators in the numerator. Recall that we need only up to 
hrst order in k since we take fc —0 limit after extracting piece. If one expands the 
denominator to linear order in k, we then need to put fc = 0 in the projection operators. 
The resulting trace using (2.16) gives 

tr [(T^Vs{p)(T^Vs{p)] = ^ , (2.20) 

where we use 

tr + g^^g^^) + 2i , (2.21) 

and pI = 0. It is clear that this contribution does not lead to a P-odd contribution 
which should be anti-symmetric in i and j. Hence, we can ignore k dependence in the 
denominator, which allows us to use the ordinary techniques dealing with pinch singularity 
in fc = 0 limit. In section 4, we will prove that this simplihcation generalizes to all order 
ladder diagrams, that is, the fc-dependences in the denominators appearing in the ladder 
diagrams do not contribute to a P-odd part of the correlation function up to hrst order 
in fc, and hence can be ignored. 

In computing P-odd fc-dependence in ti[a''Vs{p + k)a^Vs{p)] using (2.16), one can 
replace Vs{p + k) = -sa- {p + k) s / {2\p + k\) with -sa- {p + k) a j {2\p\) by the same reason 
as above, and we have 

ti[cT'T,(p + k)a’V,{p)\ ~ [cr‘(& ■ (p + J;),)cT»(a-p,)] 

~ + (2.22) 

where in the second line, we use the fact that P-odd contribution can come only from the 
last P-odd term in the a-matrix trace (2.21). Our symbol ~ cares only P-odd part linear 
in k. When p = 0, u = I we have 

where we use \p + k\ ^ |p| + {p-k)/\p\ + (P(fc^) and replace p'^p’’ —)■ (1/3)5™'*|pp since the 
angular p integration in the hnal expression (2.19) is isotropic. When p = I and z/ = 0 
we have 

-s W£"'(p' + k‘)\p\ ~ , (2.24) 

2|Pr Ap\ 
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so, summing these two possibilities gives us the P-odd part of tr [a^Vsip + k)a^Vs{p)] as 


tr [a^Vs{p + k)a^Vs{p)\ 


rsj 



(2.25) 


and from (2.19) we have 




,Pinch 


(ra) 


~ zue 


ijl 


d^p 

(dn+{p^)\ 1 ^ 


s 


(27r)4 ' 

dp^ J 3 p ^ {jp — s p 

+ </2)(p° - s|p| 

- Km 


(2.26) 


The remaining computation is a standard procedure dealing with pinch singularity 
appearing in the denominators of (2.26). The integration can be done in the complex 
plane by closing the contour in either upper or lower half plane. The leading singularity 
appears from the pole = s|p| ± iC/2 where the residue contains a factor of l/(±iC): 
this gives a leading order contribution at weak coupling limit since ( log{l/g)T. 

Once this 1/C term is identihed from the residue of the denominators in (2.26), one can 
neglect ( in the pole location ^ s|p| for all other terms as it engenders only higher 
order terms in g. This is because the p integration has its dominant support in the region 
IpI ~ T while ( ~ g‘^log{l/g)T \p\. Therefore, one can effectively replace the two 
denominators in (2.26) with 


(p° — s|p| + i(/2){p^ — s|p| 


_, 

iC/2) *C 


slpl), 


(2.27) 


which will be used frequently in the following sections. This gives us 


Pinch 

^{ra) 


c 


d^p 1 ^ / dn+{p^)'\ 


pO=s|p| 


(2.28) 


where we ignore the momentum dependence of ( for now, which is not strictly valid 
(we will be more precise in our full ladder resummation in the next section). From 
dn+{pP)/dp^ = — n+(p°)) and n+(—|p|) = 1 — n_(|p|), the integral becomes 

' dn+ (p°) \ 

/ pO—^I pI 


d^p 1 
(27r)3 3|p| 




= -/S 


d^p 


(«+(|p|)(i - ^+(|p|)) - («-(|p|)(i - ^-(IpI))) 


/S 


bTT^ 


(27r)3 3|p| 

/ OO 

d\p\ \p\ (n+(|p|)(l - n+(|p|)) - (n_(|p|)(l - n_(|p|))) = 


and we hnally have a 1-loop expression 

g«,pi„ch ^ 
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1—loop _ 






(2.29) 


/i 


bTT^ ’ 


(2.30) 



Although the overall sign of .^5 depends on the chirality, the relative negative sign com¬ 
pared to the static value of chiral magnetic conductivity (T^( 0 ) = doesn’t depend 

on chirality and is a meaningful dynamical result. Holographic computations produce the 
same negative sign between ct^(O) and .^ 5 . 


P 


k 

P + k ^P 

Figure 2: A generic ladder diagram that contributes to the leading log result. 

The above exercise shows quite a similar feature to that one hnds in the electric con¬ 
ductivity, and one can follow the lessons we have learned from the computation of electric 
conductivity. The 1/C dependence from a pair of pinching propagators Sra{p)Sar{p) sig¬ 
nals a non-analytic dependence on the coupling constant. In a multi-loop ladder diagram 
shown in Fig.2 for example, each pair of pinching propagators sharing the same momen¬ 
tum produces a factor of 1/C ~ 1 /( 5 '^ log(l/ 5 ')T) which compensates a from an extra 
gauge boson exchange, making the diagram of the same order as the 1 -loop diagram in the 
power counting of coupling constant. Hence, one needs to sum up ah multi-loop ladder 
diagrams to get a correct leading order result for Cs, which can be achieved by solving a 
Schwinger-Dyson type integral equation that we will describe in the next section. More 
elaborate power counting [52, 54] shows that the leading contribution comes from the soft 
region of gauge boson momentum Q ~ gT, so one needs a Hard Thermal Loop re-summed 
gauge boson propagator [55, 56] for the internal gauge boson exchange lines. The fermion 
momentnm stays hard ~ T, so fermion lines and ah vertices are bare ones. 

From the 1-loop result of with 1/C ~ 1 /{g‘^ log{l/g)T) dependence, the leading 

log result for C 5 from solving the Schwinger-Dyson equation might be expected to be 
~ 1 /{g"^ \og{l/g)T). However, the correct dependence turns ont to be ~ l/{g^\og{l/g)T): 
this is also the same as in the electric condnctivity. In both cases, a physics reason behind 
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this is that small angle scatterings {6 g) by transverse space-like thermal gauge boson 
excitations (whose non-zero thermal spectral density is due to Landau damping physics) 
cannot affect the charge transport phenomena much, since they deflect charged fermion 
trajectories responsible for charge transports only slightly by small angles. On the other 
hand, these small angle scatterings by ultra-soft {p ~ g^T) transverse gauge bosons is 
the dominant source for the total decay rate ( ~ 5 ^^ log(l/ 5 ')T, where the log comes 
from log(m£)/A/^) ~ log(l/ 5 f) with Ajr ~ g'^T being the non-perturbative IR cutoff 
for the transverse magnetic sector, and ~ gT is the characteristic soft scale. This 
means that the effective IR regulator for the pinch singularities that is meaningful for the 
final conductivities is not given by the total damping rate (, but is provided by larger 
angle scatterings {6 3> g) and fermion-conversion to gauge bosons, which are governed 
by g^log{l/g)T rate. In the latter, the origin of the log is completely different: it is 
from log{T/mo) ~ log(l/5')- In our diagrammatic approach of the Schwinger-Dyson 
equation, this physics manifests itself in a nice cancellation of leading log part of ( in the 
equation that we will see in the following sections, and what remains is indeed something 
of g^ \og{l/g)T coming from the rate of fermion-conversion to gauge boson. 

We end this section by recalling that the situation is quite different for color conduc¬ 
tivity where even small angle scatterings by thermal transverse gluons can change the 
color charge of charge carriers (either fermion or gauge boson) due to non-Abelian nature 
of color charges [57], so that the same rate responsible for the leading log damping rate 
also governs the color conductivity, leading to its 1 /{g"^ \og{l/g)) behavior [58, 59]. 

3 Ladder resummation of P-odd pinch singularities 
in ra-basis 

In this section, we set-up the Schwinger-Dyson equation that sums up all-loop ladder 
diagrams with leading order pinch singularities. The idea is essentially similar to the 
one in the diagrammatic computations of shear viscosity or electric conductivity [52, 53], 
except that we have to keep a finite external momentum k up to first order in k to extract 
a P-odd part (but, we can still put a; = 0 from the outset since one factor of u comes 
out from kinematics, see (2.18) and (2.19)). We choose to work in real-time Schwinger- 
Keldysh formalism in ra-basis for our convenience, rather than the Euclidean formalism 
with subsequent analytic continuation as used in some previous literature. For electrical 
conductivity, we check that they produce the same result as they should. 
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A*(p, k) 




A*(p,fc) 


r 



k 


Figure 3: Exemplar real time Feynman diagrams that can give leading pinch singularity. 

Since we are computing the vertex at the far right in any ladder diagram 

is an a-type one with one fermion leg r-type and the other a-type. Since there is no 
aa-propagator, the a-type leg should have r-type in the other end on its left. Because 
a pinch singularity can appear only from a pair of Sra and Sar, the r-type leg from the 
vertex should have r-type on the other end on its left, since having a-type on the other 
end gives the same type of fermion propagator to the one from the former, and does not 
give a pinch singularity. See Fig.3 for an exemplar ladder diagram that can give a leading 
pinch singularity. In our convention, one reads ra-types of a fermion propagator along 
the reversed direction of its momentum arrow, which can be seen in Fig.3. The reason 
why having a rr-type propagator in the diagram can give rise to a pinch singularity is the 
thermal relation (2.15), 

Srr{p) = Q “ {Sra{p) “ Sar{p)) , (3.31) 

SO that one can pick either Sra or Sar piece from Srr to have a pair of Sra and Sar that 
gives a pinch singularity. It is clear then that the rest of a ladder diagram on the left 
other than the far right vertex should have two hnal fermion legs of r-type on its right, in 
order to create a leading order pinch singularity: that is, it has to be an effective rr-type 
vertex. At 1-loop order, this was automatic since it is a bare J* vertex. What we have to 
do is to sum up all loop ladder diagrams for this effective rr-vertex that appears on the 
left side of the diagram. 

Denoting the resulting summed vertex A^{p,k) that is a 2 x 2 matrix acting on the 
spinor space, where p is the loop momentum and k is the (small) external momentum. 


11 



Sar{p) 





Figure 4: The diagrams that need to be computed to obtain the retarded response function 
to leading log order. The effective vertex on the left includes infinite number of ladder 
diagrams. 

the final obtained from two possible Feynman diagrams in Fig.4 which look 

similar to those in Fig.l except that the vertex on the left is now A*(p, k) instead of a*, 

(^jj,Pmch(^) = ('_!) J [A\p, k) {Sra{p + k)a^ Srr{p) + + k)(J^ Sar{p))] ■ 

(3.32) 

Using (3.31) picking up only pairs of Sra and Sar for a pinch singularity, and expanding 
it in u with the same manipulation that led to (2.19) gives, up to first order in ca. 

Pinch., ^ ^ _ f d^P f dn+(/) \ _ tr [A^(p, fc)P,(p + fc)n^'P,(p)] _ 

^ J (27r)4 \ dpO J ^ (pO - s\p + k\+ iCp+fc,s/2)(p° - s|p| - fCp,s/2) ’ 

(3.33) 

where we put external frequency ca = = 0 in the effective vertex A*(p, fc) and other 

places since we already have one u factor in front. Note also that the damping rate Cp,s 
depends on the on-shell momentum as well as s = ± (that is, whether it is particle or 
antiparticle) as indicated in the expression^. We will be concerned with only this object 
A*(p, k) after putting ca = 0 in the following. 

The summation of all multi-loop ladder diagrams for this effective rr-type vertex, start¬ 
ing from the bare one J* = 'ipla^'ipr can be achieved by solving the associated Schwinger- 
Dyson type equation, which is depicted in Fig.5. The “kernel” which is made of two 
internal fermion lines and one soft gauge boson (we call it photon) exchange can have 
three possible Feynman diagrams that can give a leading pinch singularity as shown in 

^The dependence on s comes via the combination sfj, in the presence of chemical potential ^ we are 
considering. See our Appendix 2 for a detailed discussion. 
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Figure 5: The Schwinger-Dyson equation for the effective vertex A*(p, k) . 



Figure 6: The three real-time diagrams with leading pinch singularity for the kernel in 
the Schwinger-Dyson equation depicted in Fig. 5. The effective vertex connects only to 
r-type endings of fermion propagators. 


Fig.6. The resulting Schwinger-Dyson equation * is written as (in the following we denote 
QED coupling constant by e instead of g) 


A\p, k) = a" + {ief j + Q)A\p + Q, k)Sra{p + Q + k)a°‘G^^^\Q) 

+ 1 ^ySarip + Q)A\p + Q, k)Srr{p + Q + k)a^G%\Q) 

+ -^^^^^rr{p + Q)A\p + Q,k)Sra{p + Q+ k)a^G^^^\Q) , 

(3.34) 

where G^^^ {a,b = r or a) are the photon propagators in Schwinger-Keldysh contour 

GP = {A<:\Q)A‘^\-Q))3 k = j<l'x e-“‘-{A‘t\x)Af(0))sK . (3.35) 

*In [53] it was shown that the Ward identity requires addition of an extra term in Fig.5 involving soft 
fermion lines. This diagram gives a sub-leading contribution to the electric conductivity and we expect 
the same for our ^ 5 . We leave the explicit computation as a future work. 
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including Hard Thermal Loop (HTL) photon self energy (that is JJ correlation functions) 
since the photon momentum Q is soft. We will work in the Coulomb gauge which separates 
longitudinal and transverse modes in a clear way. A summary of in this gauge 
including the P-odd part coming from the P-odd part of HTL photon self-energy is given 
in the Appendix 1, where we also hnd some useful sum rules for the P-odd part of their 
spectral density, which will be used importantly later. Using (3.31) and a similar thermal 
relation for photons (see the Appendix 1), 

= () + nB(y)) (Gi’jtQ) - g1“;>(Q)) = (1 + , (3.36) 

where the photon spectral density is dehned by 

P%(Q) ^ (G%'(Q) - GL“;'(<3)) , (3,37) 

and nsiq^) = l/(e^'^'’ — 1), the pinch singularity part of the integral equation becomes 

A*(p, fc) = a* + {ief J ^^[a^Sarip + Q)A-\p + Q:k)Sraip + Q+ k)a°' 

XP$(Q) (^+(P° + + nsiq^))] ■ (3.38) 

Note that the photon spectral density Pq^(Q) is hermitian in {a,/3) indices, but not 
necessarily real. In fact, the P-odd self-energy leads to a purely imaginary, anti-symmetric 
contribution to the spectral density. We refer the readers to the Appendix 1 for a detailed 
exposition. 

From the pair Sar{p + Q)Sra{p + Q + k) in (3.38) for small k limit, one can extract 
the leading pinch singularity, 

Sarip + Q)A*(p -|- Q, k)Sra{p P Q P k) 

+ q)A\p + Q, k)Vt{p + q + k) _ 

^ (pO + gO _ + q\- i(p^gt/2){p^ + q0-t\p + q + k\+ iCp+q+fc,i/2) 

Since the photon momentum Q is soft, Q \p\ ~ T, and the pinch singularity in the 
hnal equation for {k) in (3.33) necessitates the loop-momentum p to be on-shell, 

p^ = <s|p|, the only possible way to have a pinch singularity for soft q^ integration in (3.38) 
is to pick only t = s piece in the expression (3.39): it means that in a ladder diagram 
the leading pinch singular contribution comes from a particle loop or an anti-particle 
loop without any “transition” between particle and anti-particle throughout a diagram. 
Physically, it is obvious that a nearly on-shell particle (anti-particle) can not change to 
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an anti-particle (particle) with soft photon scatterings. Therefore, for a given choice of 
s in (3.33), we keep only t = s piece of (3.39) in the integral equation (3.38), and the 
solution of the resulting integral equation we also label by s: As{p,k). The more correct 
expression for (3.33) is then 


= -cn 


f ‘'‘p 1 

k dn+{p^)\ 

\ 

tr [A* (p, k)Vs{p + k)a^Vs{p)] 

/ (27r)4 ' 

\ dp^ ) 

0 

1 

0 ^ 

\P + k\ 

+ iCp+k,s/‘2){p° - s|p| 

- <p,s/2) ’ 


(3.40) 


where A* (p, k) satishes the integral equation 


KiP,k)=a^ + 


d^Q 

(27r)4 


cr^Vsip + q)Al{p + Q, k)Vs{p + q + k)a‘' 


p$(g) (n+(p° + g°)+n^(gO)) 


X 


(pO _|_ qO 


s\p + q\- fCp+q,s/ 2 )(p° + q^ -s\p + q + k\ 


^Cp+q+fc,s/2) 

(3.41) 


The rest of the paper is about solving the integral equation (3.41) in leading logarithmic 
order in the coupling constant e. 

Since our transport coefficient .^5 is obtained from the P-odd part of by 

k) = + 0(u'‘, e ), ( 3 . 42 ) 

we would like to expand (3.40) in k up to first order, focusing only on the P-odd 
structure at the same time. Since .^5 is CP-odd which shares the same quantum number 
with the (axial) chemical potential p, it can only contain odd powers in p, as seen in 
the 1-loop computation in the preceding section. In our work, we will only compute .^5 
up to linear order in p in small p limit, neglecting higher order terms of p^ and beyond. 
Therefore, we will only be interested in a linear p dependence of (3.40) and (3.41) in the 
following. 


4 An important simplification 

In solving (3.41) up to linear order in fc, and using it to compute (3.40), there are various 
sources of k dependence appearing in the equations. The problematic source is the k 
dependence in the denominators of the equations (3.40) and (3.41). For example in 
(3.40), we have 

- s\p fc| iCp+fe,s/2 ^ - s|p| -h iCp,s/2 -sp-k + i{dC,p,s/dp) • fc/2 H- (4.43) 
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giving rise to up to linear in fc, 


sp ■ k — i{d(p^s/dp) ■ k/2 


+ ■■■ (4.44) 


- s|p + fc| + iCp+fc,^/2 pO - s|p| + iCp,s/2 (pO - s|p| + iCp,s/2)2 

The second term is a double pole, and when used in (3.40) it engenders a ~ l/e de¬ 
pendence which is larger than the usual 1/C pinch singularity. The same is true for the 
fc-dependence in the denominators of the integral equation (3.41). For P-even part, this 
may be what one encounters when trying to include a hnite k in the current correlation 
functions, which seems to be related to the expected appearance of diffusion pole structure 


(^el 


UJ 




-h -\- 


UJ 




(4.45) 


since the Einstein relation gives D = (Jei/x {x is the charge susceptibility) and the k 
dependence is quadratic in electric conductivity ~ T^/e^ log(l/e). However, such a 
diffusion pole structure is not expected in the P-odd part of our interest [12] t, and it is 
natural to expect that these fc-dependences from the denominators in (3.40) and (3.41) 
do not contribute to our P-odd structure Let us show this important simplihcation 

in the following. As a consequence, one can ignore all fc’s in the denominators of the 
equations, and the only interesting k dependence comes from the projection operators in 
the numerators. 

Let us consider a generic multi-loop ladder diagram depicted in Fig.7, and let’s choose 
an arbitrary internal SraSar pair from the “side rail” that can give rise to a pinch singu¬ 
larity in small k limit. By shifting loop momentum p, the momentum that flows in Sra 
can be put to p -|- fc, then the momentum of Sar is p. The denominator of Sraip + k) 
contains a fc-linear piece as in (4.43). We would like to show that this k dependence does 
not lead to any P-odd structure 

Once we get a term like (4.44) that is linear in k from expanding the denominator 
of Sraip + fc); we should put fc = 0 iu all other parts of the diagram since it already 
saturates the linear k dependence we are looking at. These include projection operators 
in the numerators of Sraiv + k) and Sariv)'- 'k’sip + k) and Vsip) (recall that we need to 
have a same s throughout the diagram for the leading pinch singularity), as well as the 


t There could arise a diffusion pole structure in the P-odd part if one considers coupling to energy- 
momentum sector of the theory leading to “chiral magnetic energy flow”. However, when expanded in k 
it would give a term of or higher [60, 61]. Also, the coupling to energy-momentum sector is of order 
fp, and the resulting P-odd effect is of order [60, 61]. Therefore, we can ignore this possibility in our 
work. 
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Sra{P+k) 


Figure 7: Generic ladder diagram that gives rise to a leading log contribution. The ladders 
can be codified in effective vertices A,A' 


remaining parts of the diagram other than the chosen Sra{p + k)Sar{p) pair, which we 
call effective vertices: the part on the left let us call A*(p, fc) and the part on the right 
A^{p, k). See Fig.7. The value of the diagram is then proportional to (the p integral) of 


tr 


^\p)'Ps{p)A\p)rs{p) 


oc tr 


K\p){a ■ ps)A^{p){a ■ ps 


(4.46) 


where A*(p) = A®(p, fc = 0), A^{p) = A^{p, fc = 0), etc. All ij index structure comes from 
this spinor trace. 


The effective vertex A*(p) is a 2 x 2 matrix in the spinor space and since = ( 1 , cf) 
forms a basis for any 2x2 matrices, we write A*(p) = a^X^^{p). By invoking rotational 
invariance, we generally have 

Ao(p) = /o(p°, \P\)P\ AKp) = /i(p°, \p\) SI + f 2 {p°, \p\) PiP" + (4.47) 


where fi are functions only on p^ and \p\. Similarly, we have for A^{p) = (t'^XKp) with 
Kip) = /o(p°, IpDp’ , K(p) = fi(p°, |p|)5/+/2(p°, IpDpip^+Mp^, IpDK'^Pm- (4.48) 


Inserting these representations of A*(p) and A^{p) into the above (4.46), and using the 
trace formula we repeat here 


tr - g^^^g'^^ + g^K^'K 


+ 2i , 


(4.49) 


we immediately see that the last e-tensor term in (4.49) does not contribute since ps 
appears twice in (4.46), and we have 


tr 


^\p)K ■ Ps)~AKp){a ■ Ps) 


4(A*(p) ■p.)(A^(p) - Ps ), 


(4.50) 
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where we use pi = 0. Then, we have 


A*(p) -Ps = s\p\\q{p)+P^\\{p) = (s|p|/o(p°, \p\) + /i(p°, IpI) + |pr/2(p°, \P\))P\ (4.51) 


which is proportional to p*. Note that the piece involving /a drops out. The same 
conclusion is true, that is, \^{p) ■ Ps ~ and hence the result for the trace in (4.46) is 
proportional to p*pT Since it is symmetric with respect to ij, it is clear that the result can 
not contribute to the P-odd part of In summary, we have shown that k dependence 

from the denominator of any internal fermion line in leading pinch singularity limit does 
not contribute to the P-odd structure and hence we can neglect all fc’s appearing 

in the denominators, especially in our equations (3.40) and (3.41). 

Once we remove all fc’s from the denominators, we have 


Pinch 


{k) = —u 


and 


(^ 


/ dn+(p°) \ tr [Ai{p, k)Vs{p + k)a^Vs{p)] 

V dpO / ^ (P° - ^IpI + <p,s/2)(p° - s|p| - iCp,s/2) ’ 

(4.52) 


KiP,k) = (T^ + 


d^Q 

{2nY 


X 


+ q)K{p + Q> fc)'Ps(p + q + 

rfyc) (n+(p“ + 9°) + "8(9“)) 


(p» + ?» - s|p + ijI - <„+,,,/2)(p» + - s|p + q| + 


(4.53) 


The p° integration in (4.52) can be computed in leading pinch singularity limit by replacing 

1 


^ |^d(p° - s|p|), 


(4.54) 


(pO - s|p| + iCp,s/2)(p° - s|p| - <p,s/2) Cp,5 

which enforces the on-shell condition p° = s|p| on the p appearing in the integral equation 
(4.55). We will assume this on-shell condition throughout our computation in the following 
sections. Then, the integral equation becomes 

d^Q 


Ki.P^k) = a^ + e 


(27r) 


a^Vs{p + q)K{P + 9, k)Vs{p + q + k) 


a 


pfpiQ) («+(5|p| + <f) + riB{q^)) 


X 


(gO + s|p| - s\p + q\- iCp+q,s/2)(g° + s|p| - s\p + q\ + <p+q,,/2)_ 

(4.55) 


where A*(p, fc) = A® (p, fc) The reason why we can also replace A®(p -|- Q,k) in 

the integral kernel with its on-shell value A* (p -|- q, k) is that the pinch singularity in the 
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kernel of the integral equation 


- —y - 6(zP-\-Q^ — 

(pO + gO _ s|p + q\- + gO - s\p + q\+ iCp+g,s/2) Cp+q.s 

(4.56) 

will impose the on-shell condition = s\p + q\ as well. With this replacement of 

(4.56) in (4.55), we hnally have 


A*(p,fc) = a* + 


d^Q 

(27r)4 


a^Vs{p + q)K{P + 9 , k)Vs{p + q + k) 


a 


X pJW) (n+(s|p|+ 9 “)+nB( 9 “)) (2ir)i5(g" + s|p| - s|p + q|)/Cp+,,,| 

(4.57) 


which is our starting point in solving the integral equation in leading logarithmic order 
in the next section. 


5 Leading-log computation 


The integral equation (4.57) obtained in leading pinch singularity limit is a matrix equa¬ 
tion, and it is desirable to transform it into a scalar equation. In fact, we don’t need its 
full matrix structure: what we need at the end in (4.52) is the trace 

tr [A*(p, k)Vs{p + k)a^Vs{p)\ = tr ['Ps(p)A*(p, k)Vs{p + k)a^] , (5.58) 

and it is obvious that we only need the component of A* (p, k) projected onto the eigenspace 
of the projection operator Vs{p) on the left and Vs{p + k) on the right, that is, 

'P^(p)A* (p, k)Vs{p + k) . (5.59) 


Since the spinor space is two dimensional, the above projected matrix is essentially a single 
number. This fact manifests itself by the following statement; for any 2x2 matrix A, the 
projected matrix Vs{p)AVs{p + k) must be proportional to the matrix Vs{p)Vs{p + k) 
composed only of the projection operators. The proportionality constant, which contains 
the information on A, is easily found by comparing traces of the expected relation. 


Vs{p)AVs{p + k) 


ti[Vs{p)AVs{p + k)\ 

tr [Vs{p)Vs{p + k)] 
ti[Vs{p)AV,{p + k)] 
\{l+p-p + k) 


Vs{p)Vs{p + k) 
Vs{p)Vs{p + k ), 


(5.60) 
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where p = p/\p\, and tr [Ps{p)Vsip + k)] = |(1 +p p + k). Using this, it is straightfor¬ 
ward to convert onr matrix integral eqnation into a scalar eqnation, and for this pnrpose 
let us introduce three scalar functions V^p, k), k), and Fs{p', q, k) by 

'Ps{p)KiP, k)Vs{p + k) = Vl{p, k)Vs{p)Vs{p + k ), 

Vs{p)(J^^Vs{p + k) = k)Vs{p)Vs{p + k ), 

Vs{p)Vs{p + q)Vs{p + q + k)Vs{p + k) = Fs{p; q, k)Vs{p)Vs{p + k). (5.61) 


The expressions for k) and Fs{p'i q, k) can easily be found by computing the neces¬ 

sary traces involved, for example we have 


S0(p,fc) = l, KiP,k) = 


s{^ + {p + k) 


le 


ijl 


Pj{P + k)i 


(1 +PP + k) 


(5.62) 


and the expression for Fs{p]q,k) can be found in the Appendix 3. Using the fact that 
is hermitian, we also have 


Vs{p + k)a^Vs{p) = (S^(p, k)y Vs{p + k)Vs{p ). (5.63) 


The scalar function T>l{p, k) is what we would like to hnd by solving the integral equation, 
and once it is found, the hnal expression for the correlation function is given 

from (4.52) by 





tr [Vs{p)Ai{p, k)Vs{p + k)a^' _ 

(pO - s|p| + iCp,s/2)(p° - s|p| - iCp,s/2) 

Vi{p, fc)tr [Vs{p + k)a^Vs{p)] 

(pO - s|p| + iCp,s/2)(p° - s|p| - iCp,s/2) 

'k^i{P,k)m{p,k))* 4(1 + p.p + fc) 

(pO - s|p| + iCp^s/‘2){p^ - s|p| - iCp,s/‘2) 

(5,64) 

pO=s\p\ 


where we use (5.61), (5.63) in the second and third lines, and replace |(1 +p p + k) with 
unity in the last line, since p ■ p -|- fc = 1 up to negligible 0{k'^) corrections. As the last 
expression involves the combination T>*(p, k)/(p^s, let us also dehne 

(5.65) 
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Applying projection operators Vs{p) on the left and Vs{p + k) on the right to our 
integral equation (4.57), and using (5.61), (5.63) we dually get the following scalar integral 
equation to be solved for 0* (p, fc). 


Cp,s(t)\{p,k) 


KiP,k) + Sf(p,qr)0;(p + qr,fc)(S“(p + fc,qr))*F,(p;qr,fc) 

X pff}iQ) (^+(^IPI + + ^s(g°)) (27r)(5(g° + s|p| - s\p + qr|) 

Kip, k) + j fCsip, k; Q)Kip + g, fc)(27r)(5(g° + s|p| - s\p + qr|) , 

(5.66) 


with an integral kernel ICsiP, k] Q) dedned as 


)Csip,k-,Q) = Kip,q) {KiP + k,q)T p^a 0 iQ)Fsip;q 


,fc)(n+(s|p|+g°)+nB(g°)). (5.67) 


Our task is to dnd a solution for 0* (p, k) up to drst order in k which can give rise to the 
P-odd structure e^^k} in the expression (5.64) for We will be interested in 

only the leading logarithmic order in e. 

The structure of 0* (p, k) expanded up to drst order in k can be severely constrained 
by rotational invariance, 

4''Mk) = V(p) + /i'(p)fe' = X.(|p|)P‘ + a.(|p|)V + 4,(|p|)p*(P'fc) + ij^c.(|p|)e''Vfe”, 

(5.68) 

with four scalar functions Xs,cis,bs,Cs which depend only on |p|. Using the expansion of 
(S^(p, k))* up to drst order in k, 


imp, k))* = spi + ^ {5^^ - p^p^) fc' + , 


2|p 


2|p| 


(5.69) 


the expression for from (5.64) can be computed to drst order in k to dnd the 

P-odd structure proportional to as 


Pinch 

^{ra) 


{k) = —u 


d?p 

(^ 


E 


5 = ± 


3^ 


d?p 

(27r)3 


E 

5 =± 


(in+ ijP) 
dp^ 

dn+ ijP 

dp^ 


pO=s|p| 


{Kip, k))*Kip, k) 


3 = s|p| \P\ 


^ (sc,(|p|) - -x.(IpI) ) 


(5.70) 


where ~ above only cares about P-odd terms. Note that the two functions (a^, bg) do 
not contribute to our P-odd term, so we don’t need to compute them. Therefore, we will 
focus on Xs and Cg only in the following. 


21 



5.1 Computation of Xs(|p|) 


The function Xs(|p|) is something that has already been known in previous computations 
of electric conductivity, although we will see that there is an important correction to it 
linear in p that is relevant to our hnal value of .^ 5 . It satishes the integral equation (5.66) 
after putting k = 0 

Cp,.X>(|p|)p‘ = sp' + e^ j ^^K,(p, 0 ;<J)x,(|p + q|)pVq‘( 2 ir)i(?" + s|p| - s|p + q|). 

(5.71) 

Our treatment that follows for Xs(|p|) is mostly the same to that one can find in Refs.[52, 
53], and our computation for will conhrm the previous results in literature. Let us 
however present some details along which we can introduce several key elements that will 
be needed in our next treatment for c^dpl), which is new and more interesting to us. 

First it is important to observe that the integral 

that appears in the above equation is precisely equal to the contribution to the fermion 
damping rate from soft photon scatterings at full order, which contains the leading ~ 
log(l/e)T part of the total damping rate. We reconfirm this fact explicitly in Appendix 
2 including P-odd spectral density of soft photon propagator. Denoting this by Cp^’s 
(following notations in Ref. [53]), and writing the total damping rate as Cp,s = C^s + Cp^s 
where is the other remaining contribution to the damping rate arising from soft 
fermion scatterings (or more precisely, hard fermions making conversion to hard photons 
and soft fermions) which is of order log(l/e)T, the integral equation (5.71) takes a form 


Cj,sXsi\p\)p" = sp^ + e 


^s(p,0;qr) Xs{\p + q\)p + q -Xs{\p\)p" 


J (27r)4 

X {27r)5{q° + s\p\ - s\p + q\), 


(5.73) 


where we no longer have log(l/e)T explicitly in the equation, and what remains 

will be shown to be only of order ~ e^log(l/e)T. This cancellation of e^log(l/e)T 
dependence due to the identity of (5.72) with is the diagrammatic manifestation 
of the physics discussion at the end of section 2 : the relevant part of damping rate 
that is responsible for fermionic charge transport phenomena is not the total damping 
rate governed by small-angle scatterings with ultra-soft transverse photons, but the part 
arising from fermion conversions to photons with soft fermion scatterings, that is (p^. In 
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Appendix 2, we present a computation of to leading log order (see Eq.(6.204)), with 
a special care of its s/i-dependence we will need later. 

Introducing a variable = cos(p where ip is the angle between p and q, one can show 
that [53] 

+ s\p\ - s\p + q\) = ^ (1^1^ “ (^°)^) ’ 

where Q{x) is the Heaviside theta function, and 

( — loP 

Zo\q\=P-q = sq^ + -—-, (5.75) 

z|p| 

and using this, one can transform the Q integration into 

(g0)2_lq|2 ’ 

2|p| 

(5.76) 

where one needs to replace any p-g appearing in the integrand by sg°+((g°)^ —|gp)/(2|p|). 
As Q = (g°,qr) ~ eT is soft and p ~ T is hard for leading log contribution (which can be 
seen in retrospect), we expand the integrand in the integral equation (5.73) in powers of 
(g°,qr)/|p|. For example, after some algebra we have up to 0{Q^), 

x.(|p + 9l)pTA - x.(|pl)p' x.dpD + sAKIpl) + )(9°)l\;"(lp|)]) . 

(5.77) 

where we use the replacement p ■ q ^ sq^ + {{q^Y — |gp)/(2|p|) in the middle of com¬ 
putation, and x'si^) = dxs{x)/dx, etc. Similarly, we need an expansion of ICs{p,0-,q): 
with 

n+(s|p| + g°) + nsiq^) = ^ ^ 

and 

F,(p; q,0) = ^(l+p- pTg) = 1 + C>(g^), (5.79) 

what remains in ICs{p,0-,q) is the polarization-contracted photon spectral density 

ptfi{Q){^’^{p,q)T^s{p,q) , (s.so) 

where we need to expand the polarization part (S“(p, qr))* Ef (p, g) up to first order in 
q for our leading log computation. In Coulomb gauge, p'^iQ) = 0 (i = 1,2,3) and 
Pqq((5) = Pl{Q) is the longitudinal part of spectral density. The transverse part is 

= Pt{Q) + ipodd{Q)e'^^qi , (5.81) 
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where the second term is the P-odd contribution proportional to /i that arises from the 
P-odd part of current-current correlation function (or photon self-energy) in HTL limit, 
whose expression can be found in our Appendix 1. Note that it is purely imaginary, but 
anti-symmetric in i,j, so it is a hermitian matrix in i,j. For our purpose, we would only 
need its sum rules derived in the Appendix 1, Eqns. (6.155) and (6.156), 

(27r)2|g|4 + -“’ 

J-\q\ (2"^) 1^1 



up to less singular terms in small \q\ eT limit. All functions {pl, Pt, Podd) are odd 
in —)■ —g°, so we need an extra odd power of in the hnal integrand to have a 

non-vanishing integral over [—|qr|,-|-|qr|]. The only p dependence in the usual spectral 
densities Pl/t{Q) is through the Debye mass, which is m|, = e^T^/6 -|- e^p^/{2n‘^) for 
a single Weyl fermion. Since we are looking at only up to linear p dependence, we can 
safely neglect p^ corrections in m|) and use the /r = 0 results for Pl/t{Q)- After some 
algebra, we have up to 0{q) 

As(«) (s:(p, q))’ sf (p,,) = pdQ) + (i - pt(q) - , 

(5.83) 

where the last contribution from the P-odd spectral density, although it is quadratic in 
g, is presented because its power counting is something new and different from those of 
Pl/t{Q) as can be seen in (5.82), and should be checked carefully. 

Let us first estimate this contribution from the P-odd spectral density in the integral 
equation (5.73). Collecting everything presented above, the contribution from the P-odd 
spectral density to the integral in (5.73) becomes 


X 


5 1 1 

\P\ 

P‘ 

ri 

, 0 ) 2 . 

) 

- \q\ 

V 2|p 

,|2 


(27r) 


d\q\\fi\ 
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-Xsi\p\) + sq^x'si\p\) + Vs(IpI) j PoddiQ ), 


9r) 

(5.84) 


and using the sum rules derived in Appendix 1, (6.155), (6.156), it is easy to find that 
the result is at most of order ~ without any logarithmic enhancement. Note that the 
\q\ integration should have an UV cutoff ~ T since we use HTL approximation for soft 
momentum Q <^T. In any case, these are of higher order than log(l/e) of our interest. 


24 



so can be neglected. Although we find Podd(Q) does not affect the leading log equation for 
Xs(IpI), we will find shortly that it does give an important contribution to the equation 
for Cs(IpI) at leading log which is of our more interest. 


The integral equation (5.73) then takes a form at leading log order as 

d\q\\q\ dq^ ^ 

/o (27r) y_|,| W) 


Cp,.Xs(|p|) = s + e 




X ( Pl{,Q) + 1 — 
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) PT{Q)j 
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2|p|2 '^ Xs(IpI) + sgV.(lpl) + ^(g°)Vi(|p|) 


(5.85) 


The remaining computation of various integrals of spectral densities is achieved at leading 
order using well-known sum rules of Pl/t{Q) [62]. The leading log contribution will come 
from the region mD <^\q\<T, and following the notations in Ref.[53] by defining 

/-i.i WY 

we have for |g| 3> the sum rulesl- 
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Using these in (5.85), one encounters a logarithmic divergence 

r d\q\ 


= \og{T/mD) ~ log(l/e), 


(5.87) 


(5.88) 


' mo 


where the IR cutoff is ttid since the sum rule expressions used are valid only for |qr| 
niD (see our footnote), and the UV cutoff is T as we assume soft Q T throughout 
our treatment, and the modification for hard Q will dampen away the UV divergences. 
Picking up the logarithmically enhanced terms in the integral, we finally get the differential 
equation 


CpYxsHpD 


e^m|) log(l/e) 
dTT 




+ ns{\p\) 



^Xs{\p\) 

(5.89) 


^The sum rules for the case \q\ <C tod take different forms, and it can be checked that the ultra soft 
region |g| ^ mu does not give rise to logarithmic divergences. 
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which agrees with the known result in literature in the case /i = 0. Note however the 
important fj, dependence via in the differential equation for Xs(|p|)) 


M\P\) = 


^I3{\p\-sn) _|_ ]_ g/3|p| _|_ ]_ 


+ S/d/i- 




e/3|pl + 1)^ 




(5.90) 


as well as in that we compute in Appendix 2 (see Eq.(6.204)), 




dTT 
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An 


\P\ 


"fidpl) + ^ - s^j+0(p), 

( 5 . 91 ) 

which give rise to a s-independent, linear p part in Xs(|p|) addition to the usual /i- 
independent part proportional to s. Here = (e^/4)(T^ + is the asymptotic 

thermal mass of fermions. The solution when expanded in p then takes a form 


Xs(IpI) = sX(0)(|p|) + hX(i)(|p|) + , 


(5.92) 


and both X{o) and X(i) give separate contributions of the same leading order to the final 
expression for ,^5 in (5.70). In fact, this /i dependence via n^dpl) and in (5.89) (that 
is, the X(i) in the equation (5.92)) also makes a contribution to the dependence of 
the electric conductivity in leading log order, which seems to have been missed in some 
previous literature. Our analysis in the above (with full expressions for m|, and m^) 
contains all necessary elements that allow us to compute full correction to electric 
conductivity, and we present the correct computation of /i^ correction to the electric 
conductivity in Appendix 4. 


5.2 Computation of Cs(|p|) 

Let us next describe our analysis for Cs(|p|), which appears as the P-odd component of 
the 0 i(p,fc) ~ (l)l{p,0) + i{cs{\p\)/\p\)e'^''^'p^k'^ when expanded in linear k (see (5.68)), 
that satisfies our original integral equation (5.66) with a finite k. Expanding the kernel 
function ICs{p,k;Q) defined by (5.67) up to linear in k, 

}Cs{p, k- Q) = ICsip, 0; Q) + lCf\p, k-Q) + Oi^) , (5.93) 

where /Cs(p, 0; Q) is something we already use before (see Eq.(5.71) and (5.72)) to deter¬ 
mine the zeroth order solution 0* (p, 0) = Xs(|p|)p*, the part of integral equation in (5.66) 
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that is linear in k gives the integral equation for Cs(|p|), which takes the form 

IpI 

=- f -^/C,(p,0;g)^^^4^^^e*'™?Tg'fc™(27r)(5(g° + s|p|-s|p + g|) 

2\p\ J (27r)^ \p + q\ 

+ j j0^i^^P{P^^'^Q)Xs{\p + q\)p^\‘^T^)5{q^ + s\p\-s\p + q\) , (5.94) 

where in the last term we understand that we extract only the P-odd term having the 
same structure of The first term on the right arises from the fact that S* (p, k) 

given by (5.62) contains the P-odd term when expanded linear in k 


K(p,k) -— 

As before, it is important to use the fact that the integral 

/ 7W%^s{p,0;Q){27r)5{q^ + s\p\-s\p + q\), 


(5.95) 


(5.96) 


which appears in the second term on the right side is precisely equal to the contribution 
to the fermion damping rate arising from soft photon scatterings, so that one can 
transform the above integral equation into the form 

, 2 f d^Q ^ fi<^s{\p + q\) *C,(|p|) \ 


-/Cs(p,0;g) 


ics{\p + q\) 


J \ I 

X {27r)5{q^ + s\p\ - s\p + q\) 

2 f d'^Q ^(1), , . . / 


p + q\ 


+ qk^ 


iCsilpl) Jlm^lum 




+ e^y + ( 27 r)( 5 (g° + s|p|-s|p + qr|), ( 5 . 97 ) 

where C,pg that appears on the left is the damping rate contribution arising from soft 
fermion scatterings only which is of order e^log(l/e) rather than e^log(l/e). 

The computation of the first integral on the right side of (5.97) at leading log order is 
almost identical to that of the previous integral in (5.73). Expanding up to quadratic in 
g, with the replacement q^ {p ■ q)fi due to rotational invariance of the q integral, we 
get after some algebra 

Cs{\p T q\) j J^m _ *^^(|p|) Y^pYz^ 

\p + q\ \p\ 

= ( *^°!(| ~j^'' e.(|pl) + sgXm + )(9°)^g(lpl)) + o(Qt ■ (5.98) 




-Cs{\p\) + sqXm + T^idTclm e^^^p^k^ + 0{QY , (5.98) 
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where Cs(x) = Cs(x)/x. Comparing this with the previous expansion (5.77) for Xs(|p|)) 
we hnd the identical structure appearing, so that we can simply use the previous result 
of the integral in (5.73) (see Eq.(5.89)) by replacing Xs with Cs(x) = Cs{x)/x to get 
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/C,(p,0;Q) 


J (27r)4 

X {27r)5{q^ + s\p\ - s\p + q\) 
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\p\ J 2/3 


-c:(|p|))e*^Vfc™, 


(5.99) 


up to leading log order. 

What is more complicated is the evaluation of the second integral in (5.97). Let us 
hrst look at the term Xsi\p + q\)p + Q in the integrand. Dehning Xsix) = Xs{.x)/x and 
using the fact that the 5 function in the integrand imposes \p + q\ = \p\+ sq^ we have 

Xs{\p + g|)pTg* = Xs{\p\ + sg°)(K + P) ■ (5.100) 


Since what we need is the P-odd structure it is clear that the hrst term pro¬ 

portional to p* can not possibly generate such structure, and therefore it is sufficient to 
consider only the second piece proportional to g*, 

Xs{\p + q\)pTq Xs{\p\ + sq^)p , (5.101) 

in the integral of (5.97). On the other hand, since id'p{p,k]Q) is a rotationally scalar 
function linear in fc, rotational invariance dictates that it can only have three possible 
structures 

k- Q) = (a^p' + h^pq^ + cPe^'^^fTqp , (5.102) 

where (ai^\ ci^^) are some coefficient functions that depend only on {p^,\p\,q^,\q\-,p ■ 

q). Combining these two facts, and considering rotational invariance of q integration, one 
can easily hnd that the only way to have the resulting P-odd structure from the 

second integral in (5.97) is via the third term in (5.102), that is, we only need to hnd the 
part of lcP{p,k-,Q) that is proportional to = ^imnpiqmj^n^ This simplihes 

our computation by a great amount. 

Since (5.101) is already linear in Q, for a leading log contribution we only need to 
expand k!'s\p, k; Q) up to linear in Q which is already saturated by the wanted structure 
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This in turn implies that one can neglect sq^ correction in (5.101) to have 

Xs{\p + ^ Xs(|p|)g* = ’ (5.103) 

in the integral of (5.97). Given the expressions for S*(p, k) and Fs{p', q, k) in (5.62) and 
Appendix 3, as well as the photon spectral density given in (5.81), 

Pm = Pl{Q) , pf^iQ) = Pt{Q) (^Sij - + ipoddiQY^^qi , (5.104) 

it is straightforward to hnd after some amount of algebra that 
^f\pY\Q) = {n+Yp\ + q^) + UBiq^)) 

(5.105) 


up to linear in Q, which will contribute to the leading log result of the integral in (5.97). 
Note that we have a non-negligible contribution from the P-odd part of the spectral density 
Podd(Q): from its sum rule given in (5.82) one can easily see that this term engenders a 
leading log contribution to the integral. When combining lcY{p,k',Q) with (5.103) in 
the integral of (5.97), one has to replace qr”^qr* with 


^ = h”* (|q|2 - (p . qf) = 1 (|q|2 - (q«f} , 


(5.106) 


where qr is the perpendicular component of q to p, and we use (5.75) in the last equality. 
This comes from the rotational invariance of q integral around p axis. Collecting all these 
and following the same steps as in the leading log computation of y^dpl) before, we dually 
have the second integral of (5.97) to be given by at leading log order 
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where we use the sum rules (5.87), (6.155), (6.156), and interestingly it turns out that 
the contributions from the P-even spectral densities pl/t cancel with each other exactly. 
We don’t have a good understanding whether this has to be the case by some symmetry 
reason or it is just by accident. Therefore, the only contribution to the second integral 
in (5.97) at leading log order (that is, e^log(l/e)) comes from the P-odd part of the soft 
(HTL) photon spectral density Podd(Q)- Note that for this contribution, we have equal 
logarithmic contributions from both e^T \q\ ^ uid and T 3> \q\ 3> that add up 
together in the hnal result 

From the integral equation (5.97) with (5.99) and (5.107), we hnally obtain the sought- 
for second order differential equation for Cs(|p|) as 


Cp,sC,(|p|) 


1 log(l/e)/ix,(|p|) 

2 (47r)27r2 (3\p\ 


e^m\) log(l/e) 
An 





Cs(|p|) \ 

\P\ ) 



(5.108) 


where the hrst line is an inhomogeneous source, especially the second term is in terms 
of Xs(|p|) that should be obtained by solving the differential equation (5.89). We would 
need the expansion of Cs(|p|) up to hrst order in chemical potential p, 

Cs(IpI) = C( 0 )(|p|) + sp C(i)(|p|) + 0{p‘^) , (5.109) 

which can be found by solving the above differential equation order by order in p. We 
reemphasize that there are linear p dependences coming from and ns(|p|) in the above 
differential equation, which should not be missed to get a correct leading log answer. After 
hnding Xs(|p|) and Cs(|p|) from the above given differential equations, we compute our 
transport coefficient .^5 by (5.70) with 


5.3 Numerical evaluation 

As a hrst step to compute the explicit value of .^ 5 , we solve numerically the equations 
for (Xs(IpI), Cs(IpI)) in (5.89) and (5.108). In order to do so we dehne T^dpl/T) = 
amf) log(l/e)cs(|p|)/T and "h^dpl/T) = Q;m|, log(l/e)xs(|p|)/T where a = e^/ 47 r. Dehn- 

^Strictly speaking, QED does not possess ultra soft magnetic cutoff ~ e^T. Since we have in mind 
the generalization to non-abelian QCD discussed in section 5, we simply assume this at this point. 
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ing y = |p|/T, the equation (5.108) can be rewritten as 
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In addition we have the equation for the even vertex, Xs(|p|) (5.89) 
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(5.110) 
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sech2(l//2)$;(?/). 


(5.111) 


In the above, we expand y dependence from and n^dpl) up to linear in y. We then 
expand the solution to hrst order in y as 


■^s = '^A + s^^B + 0{y^), 

$, = S$A+|$S + C^(/i'), 
from which we have a coupled set of differential equations 


(5.112) 

(5.113) 
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(5.114) 


We solve the above equations by iterative method with vanishing boundary conditions 
at the IR {y = 0) and the UV {y = cx)). The last step is then to obtain the expression for 
the transport coefficient ^5 as an integral of the above quantities. After computing the 
sum over s = ± and performing angular integrations in (5.70), we obtain the result for 
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the retarded propagator as an integral of ^A/siy) and ^A/siy)) 




Tne^ log(l/e) 
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2 J cosh^(|//2) J ■ 

(5.115) 


From the identihcation we hnd 


3.006 /i 
log(l/e) T 


(5.116) 


6 Discussion 


It is easy to generalize the above to the case of Np species of Dirac fermions in an SU{N^) 
gange theory (but still the current and magnetic held are with respect to the global U{1) 
havor symmetry). The chemical potential appearing in the integral equation is simply the 
axial chemical potential ^a- The Debye mass and asymptotic thermal mass are changed 
to 

o 2 t ^2 2rp2 7^2 _ 1 

y AT , AT \ 2 y c 
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i2N, + NF) 


mf = 


2Nr 


(6.117) 


The soft-fermion contribution to the hard photon damping rate is given by 


Asf _ y ^''c 


1 m} \og{l/g) 


Att 2Nr 
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(^i?(|p|) +^-s( 0 )) 


(6.118) 


In the integral equations (5.89) and (5.108), the e^log(l/e) in the kernel part should be 
replaced by 

- 1 

y^iog(i/y)^^. (6.119) 

In the second source term in the hrst line of (5.108) which is proportional to e^log(l/e), 
this e^log(l/e) should be replaced by 


^og{l/g)NF 


N!-l 

2N, 


( 6 . 120 ) 


This is because one factor of coming from the spectral density Podd(Q) is replaced by 
g'^Np, while the other coming from fermion-gluon couplings in the kernel is replaced 
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by — l)/(2iVc). Finally, the expression for {k) from (5.70), or equivalently 

our ,^5 from Xs(|p|) and Cs(|p|) must be multiplied by the fermion degeneracy 


2N,xJ2QI, 

F 

where Qp are charges of F flavors in units of e (for (u,d)-quarks, it is Qu = 
Qd = —1/3). For 2-flavor massless QCD = 3) with = 2/3 and Qd = — 
result for ,^5 then becomes 

tQCD _ 2.003 /i 

“ ^7Mog(l/^7) T • 


( 6 . 121 ) 


2/3 and 
1/3, the 

( 6 . 122 ) 


As mentioned at the end of section 2 , the color conductivity cTc for non-Abelian gauge 
theory that appears in the low energy effective theory at the scale Q < 


, 


(6.123) 


where a denotes adjoint color charge and is the thermal noise via fluctuation-dissipation 
relation to cTc, is governed by scatterings with ultra-soft transverse thermal gluons of mo¬ 
menta Q gT with the rate ~ g‘^\og{l/g), leading to <Jc ^ g^ ■ 1 /( 5 '^ log(l/ 5 f)) ~ 
1/log(l/ 5 ').^ To get a correct leading log result for this, one also needs a similar dia¬ 
grammatic resummation with essentially the same technique in our computation, except 
that charge carriers include gluons as well as quarks, and there is now no longer a pre¬ 
cise cancellation of 5 '^log(l/ 5 ') terms in the integral equation: this is due to the absence 
of U(l) Ward identity (replaced by non-Abelian version of Slavnov-Taylor identity) that 
ensures the cancellation of 5 '^log(l/ 5 ') terms [53]. The diagrammatic resummation (done 
in Ref. [63]) for this is therefore somewhat simpler than that for the electric conductivity, 
since one does not need to go to the next order of g‘^\og{l/g). In fact, one gets an alge¬ 
braic equation to solve rather than a differential equation. The same resummation can 
also be achieved in the language of Bodecker’s approach [58] as well as in kinetic theory 
[59]. In the presence of axial charge fiA breaking CP symmetry. Ref. [64] recently obtained 
via Bodecker’s approach a CP-odd contribution to the color current 

= acE'^ + e + , (6.124) 

47r^ 

which is a colored analogue of chiral magnetic effect consistent with the U{1)aSU{“ iYc 
triangle anomaly. Since triangle anomaly is topological, this contribution should be sat¬ 
urated at 1 -loop diagrammatically without a need for resummation of ladder diagrams. 


^Note that we now put an extra 
literature. 


in the definition of color conductivities to follow the convention in 
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From our computation in the text, the quantity that needs a ladder resummation and is 
sensitive to the same ultra soft scale dynamics of g‘^\og{l/g) rate that the color conduc¬ 
tivity is also sensitive to, appears when one goes to the next order in derivative 






f^ANpg^ a I 


It is clear that a colored analogue of our ,^ 5 , will be of order 


^5 ~ ■ l/(/log(l/^))(/^A/T) ~ (l/log(l/^))(/iA/T), 


(6.125) 


(6.126) 


due to the absence of precise cancellation of 5 '^log(l/ 5 ') terms in the integral equations. 
The computation of ^5 at leading log order is doable, following the same steps we present 
in our work keeping only g"^ log(l/ 5 ') terms in the integral equations (note that it receives 
contributions only from quarks, not from gluons). The leading order fermion damping 
rate from soft gluon scatterings is 


= 

^p,s 


N^-lgHog{l/g)T 


2iVe 27r 

and the solution of the integral equations which become algebraic is 

Xs(|p|) = ^ ^ g^log{l/g)T = /V^j2 ( If*!) = “ 2 ^ ’ 

which gives our result for ,^ 5 , 

g'^Np f (fp drtpijP) 
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(27r)^ dp^ 
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E fsc,(|p|) - E.(IpI) 

p^=s\p\ \ 


3n Nc log{l/g) T ' 


(6.127) 


(6.128) 


(6.129) 


The same resummation should also be achievable in the Bodecker’s approach presented 
in Ref. [64] by going to the next order in time derivatives. 

It is utmost important to implement the correct value of chiral magnetic current in 
the presence of time-varying magnetic held in any realistic simulation of chiral magnetic 
effect (or any other anomaly induced transport phenomena) in heavy-ion collisions. Our 
result should be an important step at weak coupling picture toward taking into account 
time-varying nature of the magnetic held in heavy-ion collisions, and will be instrumental 
in the quantitative studies of the chiral anomaly induced phenomena in the experiments 
at RHIC and LHC. 
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Appendix 1: Sum rules for the P-odd part of HTL 
photon spectral density 

Let us start from the thermal relation 

G" («) = (i + nB(y)) (G™(«) - G-(,)) , (6.130) 

where = (A“(a;)A^(0))5if (a, 6 = r,a) are correlation functions in the Schwinger- 

Keldysh path integral, and 

G%{x) = {A<;,{x)A:{0))sk = Glli-x), (6.131) 

where we use translational invariance of the system. Since what we encounter in writing 
down our integral equations in the main text is the combination (G™(g) — G“0(g)), let 
us naturally dehne the photon spectral density (including possible P-odd contributions in 
general) 

= G~(«) - G-(,). (6,132) 

We will show that pf^{q) is in general a hermitian matrix in terms of pu indices. For diag¬ 
onal components that come from the usual P-even contributions, is therefore 

real. For P-odd contribution which turns out to be anti-symmetric in spatial ij indices 
(there is no P-odd contribution to time-like component, at least up to linear order in p), 
we thus have purely imaginary. 

To show that p^{q) is a hermitian matrix, recall that the usual retarded propagator 
is dehned as 

G^^{x) = -ie{x^){[A^{x),A^{{))]) = -fG™(a;). (6.133) 

It is not difficult to show, using the hermiticity of A^, that Gj^^(x) is real valued, and this 
is what it should be since the retarded propagator gives the response of the system in real 
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time which must be real valued. Therefore, in Fourier space, one has 


G* (-9) = (G« (,))• , (6.134) 

which in turn gives 

G~(-9) = - (G~(9))* . (6.135) 

On the other hand, from (6.131) we have 

G;;(9) = G';(-9) = - (G~(9))‘ , (6.136) 

where we use (6.135) in the last equality. Therefore 

PlKl) = G~ (9) - G-(9) = G~(9) + (G™(9))* . (6.137) 

which proves that pf^{q) is indeed a hermitian matrix. 

The P-odd part of the retarded current-current correlation functions, that is the re¬ 
tarded photon self energy, in Hard Thermal Loop (HTL) limit has been recently computed 
in literature [65, 66]. We will work in the Coulomb gauge where = 0 (i = 1, 2, 3). The 
P-odd contribution appears only in the spatial transverse part of the correlation functions, 
so we will discuss only the spatial transverse part of current correlation functions in the 
following. In matrix notation, the spatial part of the HTL resummed photon propagator 
is 

(G’'“(9))-' = (g;;,( 9))-' - ffl«(9), (6.138) 

where G(o)(q') is the bare propagator which is given by 

The HTL self-energy S^(g) including P-odd contribution is given by 

Eg(g) = nr(g)P’5^(g)+modd(g)e*V, (6-140) 

where 
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with m|) = e^(T^/6 + /i^/(27r^)). From this while keeping terms only up to linear in p,, 
we have 


/^ra 

<^ij 


( 9 ) = 




+ 


nodd(<?) 


(g0)2 + |qr|2 - IlT{q) (-(g°)^ + |gP - nr(g)) 




(6.142) 


where the hrst term is the usual P-even HTL photon propagator, and the second term is 
the new P-odd contribution. The HTL photon spectral density pf^{q) is then given by 


= PT{q)P^j{q) + *Podd(g)e*^'g', (6.143) 


with 

priq) = 2Im 


-(g0)2 + Iql^-Uriq) 


Podd(g) = 2Im 


nodd(Q') 


(-(g0)2 + |qr|2 - nT(g))^ 

(6.144) 

It is easy to see that Podd(Q') is an odd function in and what we need in the main text 
is the value of the integral 


r\Q\ 

^(gT"-Vodd(g), n = 0,l,2,---. (6.145) 

J—\q\ 

One can compute them using the well-known sum-rule techniques exploring analytic prop¬ 
erty of the function Aodd(Q') dehned by 


'^odd(*?) 


nodd(Q') 

(-(g0)2 |g|2 _ TlT{q)y ' 


(6.146) 


We briefly sketch the procedure and present the results in two different regimes |g| mn 
and \q\ 3> niD- 


The starting point is the fact that Aodd(Q') in the complex plane is analytic in the 
upper half plane due to the causal nature of a retarded function. Thus, the integral 


dq^ 


l_^ 271 q^ — u + ie 


Aodd(<?) 0 , 


(6.147) 


vanishes for any real number uj by closing the contour with the upper hemi-circle at 
inhnity (and Aodd(Q') —t 0 sufficiently fast as |g°| —)■ cx)). From l/{q^ — u + ie) = Vl/{q^ — 
u) — i7iS{q^ — u) where V is the principal integration, we have 

r°° da° 1 i 

V / ^^^Aodd(g) - X Aodd(a;, |g|) = 0 . (6.148) 

J—00 Q OJ 2 
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Considering the imaginary part of the above, we obtain one of the Kramers-Kronig dis¬ 
persion relations for a retarded fnnction (the real part gives the other dispersion relation), 


'■oo TO 


V 


dq^ 


-Podd{q) = Re [Aodd(c<;, \q\ 


1-00 ‘271 q^ -u 

Setting a; = 0 and using Aodd(0, |g|) = — e^/i/(47r^|qr|"^), one gets a sum rule 


(6.149) 


r ^ 1 t ^ - 




47r2|qr|^ 


(6.150) 


Other sum rules are obtained from (6.149) by expanding both sides in ca —)■ cx). The 
left-hand side becomes 


-E 

n=0 


'■oo TO 


n+l 


^(g°)>odd(g), 


(6.151) 


while the right-hand side when expanded in large u is — e^/i/(127r^a;^) -|- C>(l/a;®), which 
gives other two sum rules. 


^^dq^ 0 


Tj-g Podd(g) = 0, / — (g ) Podd(g) = 


27r 


27r 


e^jj, 

127r2 


(6.152) 


The sum rules (6.150) and (6.152) are not precisely defined in (6.145), as the 
integration range for is [—|g|,-f |g|], not [—cx,-|-cx)]. The imaginary part of A(g) 
(that is, Podd(g)) consists of two distinct parts: one part coming from a branch cut just 
below the real line along the interval G [—|qr|,-|-|qr|] from the logarithms in (6.141) 
(originated from the Landau damping), and the other part is from the two poles icuo 
satisfying —(a;o)^ +Ig'p —n'r(n;o, |g|) = 0 corresponding to the time-like transverse photons 
in the medium. The former has a continuous support in the interval G [— |g|,-l-|g|], 
and hence contributes to while the latter pole contributions sit outside the interval, 
uq > |g|, so do not contribute to Therefore, the only difference between the sum 

rule values in (6.150), (6.152) and the is simply the latter pole contributions which 
we can compute. 

Near the pole location q^ uq — ie, we have the expansion 


(g°)2” ^2Aodd(g) ~ ^ ^+ {regular} , (6.153) 

(g^ — a; -I- ze)^ q'^ — Uq + le 

where 

Q ^ 2(c^o)^”'~^nodd(a;o) / nodd('^o) i‘2n - 1) _ 2 ny(a;o) \ ^ 

(2c<;o + n^(i^o))^ vnodd(i^o) 1^0 2uJo + n}T(a;o) / 
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with = dY[T{,(f )\q\)!dq^, etc. The hrst double pole does not contribute to the 

imaginary part in e —)■ 0 limit while the second part contributes to Podd^q) as 

—tiB 6{q^ — cuo) X 2 = —27rB d(q^ — coq) (the factor 2 comes from having two poles icuo), 
leading to the difference between the sum rules values in (6.150), (6.152) and the 
being given by —B, that is, is obtained by adding B to the sum rule values in (6.150), 
(6.152). 

In the case |qr| niD, the pole location is cjq ~ a/1/3 m^(l + (9/5)|qrp/m|)) + ■ ■ ■ 
and an explicit computation of B gives the values of in this regime as 


Todd ^ 
Jq ~ 


odd 


e^u 3e^u ,, , ^, 


r 

Todd 

Jo 


3e^/i|gp 

33e^p|g|^ 
' 70071 ^ 711 % 


+ 0{\q\Vm%), 


+ 0{\q\yml) 


On the other hand, in the regime |g| 3> mo, the poles are located in cuo 
(l/4)m|,/|g| + • • •, and we have the results in this regime |qr| 3> niD as 


(6.155) 
s |g| + 


j; 


jodd ^ 


j: 


odd _ eV 

STT^Igl" 
e^/i 

STT^Igl" 

- - 


1 + log 


3 + log 


m 


D 


m 


D 


m 


D 


(6.156) 


The above results (6.155) and (6.156) will be used in the main text in section 3. 


Appendix 2: Hard fermion damping rate 

In this appendix, we compute the damping rate of hard fermion including possible depen¬ 
dence on the chemical potential p. Our primary objective is two-fold: we hrst would like 
to conhrm that the integral that we have in (5.72) is indeed precisely equal to the damping 
rate induced by soft photon scatterings at full order in e and p, which was instrumental 
in rewriting the integral equation to take the form (5.73) that contains only e^log(l/e) 
terms. Our second objective is to hnd a linear sp dependence of the (p ^, that is, in the 
damping rate induced by soft fermion scatterings (or equivalently, fermion conversion- 
to-photon processes). This sp dependence in (p^^ is important in hnding the correct p 
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dependence in the solution of differential equation (5.89) for Xs(|p|), which is crucial to 
get the correct result for ,^5 as well as [j? correction to usual electric conductivity. 


gS''{Q) 




Figure 8 : Real-time Feynman diagrams for the retarded fermion self energy. 

Let us start with the self energy resummed ra-propagator (which is equal to i times 
of retarded propagator) 

^(p) = y5(p) + ss(p)E™(p)y„Vp) + ■ ■ ■ = . ( 6 . 157 ) 

which gives 

{S^%p))-^ = - E™(p), (6.158) 

where 

«(»)<?) = E . (6-159) 

Q ir \t'\ 

is the bare ra-propagator, and the self-energy Y7°‘{p) (which is a 2 x 2 matrix in the spinor 
space) at naive lowest order in coupling is given by two Feynman diagrams in Fig .8 in 
real-time formalism with the expression 

E''“(p) = (-ie)V'’ j [G;yQ)S~(p + Q) + GZ(Q)S"(p + Q)] , (6.160) 

with the photon propagators G“^((5) {a, b = r, a) (see the Appendix 1 for our notational 
conventions). In the above expression, we haven’t specihed whether the propagators 
appearing in the loop are bare or HTL resummed ones, since depending on the situations 
we can consider different approximations for them to get the right leading order quantities. 
For example, if the external momentum p is soft and one is interested in the Hard Thermal 
Loop (HTL) approximation, it is enough to consider hard loop momentum Q and the both 
propagators in the loop are the bare ones. On the other hand, in the case of damping rate 
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with a hard momentum p, which is proportional to the imaginary part of the self-energy, 
the leading contribution comes from when one of the two loop propagators carries soft 
momentum (that is, either Q or p + Q), and the soft propagator must then be the HTL 
resummed propagator while the other hard propagator is the bare one. 

Rotational invariance dictates the self energy to take a form 


E™(p) = A{p^, \p\) + B{p^, \p\)p-a = iY^ Sf (p)P,(p), (6.161) 


5 = ± 


where 


Sf(p) = -i (dl(p°, IpI) + sR(p°, IpD) . 
From this and (6.158) we have 

i 




pO - s|p| + S^(p) 


Vsip) 


(6.162) 


(6.163) 


In deriving above, we use the following properties of the projection operators to hnd the 
inverse of 5'’’“(p), 


V+ + V. = l, Vl = V±, V+V.=V.V+ = 0. (6.164) 


By comparing the above expression for S^°'{p) with the one in (2.11), we see that the 
damping rate is given by the imaginary part of retarded self energy Sf (p) at on-shell 
momentum p° = s|p|. 


Cp,s 


2Im [Sf(p) 


pO=s|p| 


(6.165) 


In the following, we will hence only concern about the imaginary part of Sf(p). From 
(6.161) and tr['P±] = 1, we have 


Kip) = {-t)tT[Vs{p)^^ 


= ie^ti Vs{p)<J^ 


d^Q 


[G':s(Q)S'^‘‘(p + Q) + + Q)] a‘ 


(6.166) 


which will be the starting point of our computation. 

For a soft p, if one uses the HTL approximation to the retarded self energy S^(p), the 
result is the HTL fermion propagator. For a hard p, the HTL self energy is sub-leading in 
so can be negligible, and moreover its imaginary part at on-shell momentum p° = s|p|, 
which would give a damping rate that could regularize pinch singularities, vanishes due 
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to kinematic constraints. The leading contribution to the imaginary part of at 

on-shell momentum arises when either Q or p + Q is soft, so that the corresponding 
propagator in the loop is the HTL resummed one. Calling the case of soft Q the soft- 
photon contribution, and the other case of soft p + Q the soft-fermion contribution, 
the total damping rate is the sum of the two, = C^s + Cp,s- 

Let us discuss first. Since the fermion propagator is the bare one, we have after 
putting on-shell momentum p^ = s|p| 


S™(p + Q) = 


^ s|p| +q^ - t\p + q\+ie 


-Vtip + q) 






P 


s|p| +q^ -t\p + q\ 


ttS {s\p + q° - t\p + q\) 


'Pt{p + q), 

(6.167) 


where P denotes principal value. Similarly, 

S^^{p + g) = Q - n+{s\p\ + ^(27r)5(s|p| + - t\p + q\)Vt{p + q) . (6.168) 

Looking at the structure of (6.166), we have a spinor trace appearing 

tr \Ps{p)(T^Vt{p + g)a“] = i7f“(p, q ), (6.169) 

which is a hermitian matrix in terms of a/? indices (this can be shown easily using her- 
mitian nature of and V±)- Since 

GUQ) = () + Ml”)) KW) - GSp(Q)) ^ () + Ml”)) pffi.Q ). (6.170) 

is also a hermitian matrix with {G'^p{Q))* = —(Q) as shown in the Appendix 1, we 
see that Hg^{p, q)G^J^{Q) is a real number. Therefore, one sees that the imaginary part 
of Sf (p) given in (6.166) arises only from the second (5-function term in (6.167) when 
used in the hrst term of (6.166). Similarly, from 

(//f”(p,g)G:yQ))' = -M.t{p,q)Gl%Q ), (6.171) 

and 

= {PnP,<l)G''MQ) - Pr(p,q)GZ{Q)) , (6.172) 

we see that the real part of P[g^{p^q)G^p{Q) is equal to — (l/2)i7fi“(p, qr)pP^(g). With 
(6.168) the imaginary part of S^(p) from the second term in (6.166) only comes from the 
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real part of Hg^{p, q)G^i^{Q), and therefore we can effectively replace appearing 

in (6.166) with —(l/2)pj(^((5) for the purpose of damping rate computation. One then 
observes that the pieces in (6.167) and (6.168) that contribute to the damping rate are all 
proportional to the (5-function (5(s|p| -|- ~Ap^Q\) which has a non-zero support only 

for t = s since Q is assumed to be soft while p is hard. After collecting all these pieces 
contributing to the imaginary part of Tif{p), we dually have after some algebra, 

CZ = 2Im[Sf(p)] 

= + n+{s\p\ + g°)) H^ip, g)p5(Q)(27r)(5(s|p| + - s\p + q \), 

(6.173) 

where 

H^^{p,q),,=tT[Vs{p)a^Vs{p + q)a<^] = Sf (p, g) (S“(p, g))* ^ (l + p ■ ^) , 

(6.174) 

using the notations introduced in (5.61). From the fact that Fs{p;q,0) introduced in 
(5.61) is equal to ^ ^1-|-p ■ p-|-, and recalling our definition of kernel function in 
(5.67) 

ICs{p,0;Q) = sf(p,qr) Q)T P$iiQ)Fs{p; q,0){n+{s\p\ +g°) +nB(g°)) , (6.175) 

we see that (Z is indeed equal to 

CZ ^ ^)(27r)<5(s|p| + g° - s|p + qr|), (6.176) 

which is precisely what appears in (5.72) and in the integral equation, which is crucial to 
have (5.73). 

Although we don’t need the value of (Z i'^ i'his work, it is easy to compute it from 
the above expression. From (5.76) 

2|p| 

(6.177) 

and the small Q expansion of /Cs(p, 0; Q), 

/C,(p,0;g) ^ ^ (^Pl{Q) + Pt{Q) (l - , (6.178) 
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;(27r)(5(s|p| +g° 




(27r) i_u|(27r) 


1 + ?^ 


I p-q^sq^ 



where Pl/t{Q) are P-even longitudinal and transverse photon spectral densities dehned 
by 


pS(Q) = pl(Q) , = pt{Q) (% - 




QiQj 


+ iPoddiQ)^^^^ 


qi 


we have 


c 


sp 

p,s 


d\q\\q\ 

(27r) 


/_|g| (27r) /3g0 


Pl(Q) + Pr(Q) ( 1 




0'i2 


(6.179) 


(6.180) 


The rest is to use the sum rules for Pl/t{Q) that can be derived by the same way we 
derive the sum rules for the P-odd part in Appendix 1 [62]. The leading log arises from 
the momentum region |qr| mo and only from the transverse part for which we have 


‘‘‘1° 1 /n\ 1 

Pt{Q) = T-r‘ 


J-\q\ (27r) g' 
and this gives 


+ 0 ^ 


m 


D 




-kl 


dq^ 

W'. 


Pt{Q) — 


3|gP 


+ 0 
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D 


pmo 


CZ ^ 


d\q\\q\ 


'>p,s 


Id 


log(l/e)T 

2^r 


(6.181) 


(6.182) 


/e^T (27r) |g|2 

where we put an IR cutoff of order Air ~ e^T. Strictly speaking, the e^T (or g^T for 
non-abelian theory) magnetic conhnement scale exists only for non-abelian theory, while 
an abelian QED which becomes free at Q ruf does not possess any IR cutoff. In this 
case, the damping rate is not a useful concept [67], and the effective IR cutoff is 
provided by the time-scale one is looking at, so the hard fermions decay in time t as [67] 




{rriDt) 


- V-t 


AlR~l/t 


(6.183) 


Since we are ultimately interested in QCD (see our discussion in section 6), we don’t 
worry about this any more. Another aspect is that in realistic situations, the free nature 
of QED at Q eT means that this scale is never thermalized anyway. Since the damping 
rate arises from the scattering of fermion with thermally excited soft photons in this scale 
(recall nsiq^) ~ l/g° term in the above), we wouldn’t have these contributions in realistic 
situations in any case. This also justihes our use of Ajr ~ e^T in the above. 

Let us next compute the soft-fermion contribution to the damping rate, with our 
main objective being to hnd a linear sp-dependence. Since p + Q is soft, it is convenient 
to shift the loop momentum Q ^ Q — p to have 


Sf (p) = ie^ tr 




d^Q 

(27r)4 


[GZiQ - P)SZQ) + GZiQ - 


(6.184) 
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where now Q is soft, and we need to use HTL resummed fermion propagators while 
the bare propagators are used for photon propagators. The HTL resummed fermion 
ra-propagator is written as 


S™(Q) = 


,0 _ i|q| + 


niq) , 


(6.185) 


where is the HTL self-energy. An explicit computation gives (in fact, one uses 

the same expression (6.184) with now both Q and p + Q being hard) 

/ / „o 




(Q) = 


m 


f 


4|g| 


2t + 


1-t^] log 
1^1 



1^ 

+ ie 

1 

O 

\q\ 

+ ie 


where 


mi = -{T^ + ^ 


(6.186) 


(6.187) 


is the asymptotic thermal mass of fermions for a single Weyl fermion system. Since 


p dependence is only quadratic for we can ignore this dependence in the 


HTL self-energy to use p = 0 result of Sf’^'^'"((5). This means that and 5'“'’ = 

— (S"’“((5))^ can be replaced by their p = 0 values up to linear in p of our interest. On 
the other hand, the rr-propagator which is given by 


S"{Q) = ( ^ - n+(9”) I («’■“(« - S"(«) . 


( 6 . 188 ) 


does contain a linear p dependence via its statistical factor in front, n+(g°). We will 
indeed hnd shortly that this will be the (only) source of the hnal s/i-dependence of 

Since the HTL resummed ra-propagator (6.185) is analytic in the upper complex 
plane, one can introduce real spectral densities p±{Q) by 


duj ptiuj,q) /■“ 


du pt{uj,q) 

gO-t|qr| + Sf’^^^(Q) J{2ti) - uj + ie ^ J{2n) - u 

or equivalently 

PtiQ) = -2Im 
in terms of which we have 


Introducing 


gO-t|g| + S™(g) 

= Q -n+(g°)^ ^p,{Q)Vt{q). 
L4“(p,g)^tr[iP,(p)a^iP,(g)a“] , 


(6.189) 

(6.190) 

(6.191) 

(6.192) 
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which is a hermitian matrix, a similar discussion to that we have above leads us to replace 


S™(0) ^ ^Pt(Q)Vt(q ), 


(6.193) 


for computing the imaginary part of in (6.184). 

On the other hand, the hard photon propagators in (6.184) are bare ones. In the 
Coulomb gauge we have 


GSo(p) = 


\P\ 


Gl^ip) = 




— (p° + iey + IpP ’ 


(6.194) 


where {p) = Sij — PiPj is the transverse projection operator, from which we have the 
bare photon spectral density as 


doo(p) = Pfip) = (2vr)P’^(p)sign(p°)(5 ((p°)2 - |p|2) , (6.195) 

with G^Ji^{p) = (1/2 + ns(p°)) p)(^(p). As before, for the imaginary part of Sf(p) in 
(6.184), we can replace 

- p) ^-^p'a^iQ - p) ■ (6.196) 

Collecting all these elements, the expression for becomes 

C!* = ^ P'^ + ^+(^°)) di(Q) 

X sign(g° — s|p|)(27r)(5((s|p| — g°)^ — Ip — qrp). (6.197) 


Since Q is soft while p is hard, we have sign(g° — s|p|) = —s and 

6{{s\p\ - qy - Ip - = 2 (|p| “ ^\p ~ ^1) ’ (6-198) 

and using the identity 

- -s|p|) + n+(g°) = (-s) (ns(|p| - sg°) + n_^(-sg°)) , (6.199) 


we have 


C = 

t=i 


X 


dy 


LitiP, q)Pij{fl - p) (^s(|p| - y) + n-s{-sqy pt{Q) 


2(|p| - sgO 


(27r)4 

(27r)(j(s|p| -q^ -s\p~q\) 


( 6 . 200 ) 
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It is straightforward to compute the leading log part of the above integral by expanding 
the integrand in powers of Q/T or Q/|p|, the same procedure we use several times before. 
From Lil{p,q)Pl^{q - p) ^ 1 - st^ + 0(Q) = 1 - + 0(Q) and ns(|p| - sg°) + 

n-s{—sq^) ~ ubUpI) + ri_s(0) + 0{Q), the leading log comes from the expression 


^ 91^ (^s(IpI) + n_,(0)) 


d\Q\\Q\ 




2|p| 

and using the sum rulesll 


271 J_|q| 271 


( 6 - 201 ) 


, m'l [ 4|qp \ 

Jo = log ^ - 1 , '^1=^,1 
4|q|2 I ^ m2 / ’ 4|q 


^'log4f-3|. 


m 


f 


( 6 . 202 ) 


where 


we hnally have 




dq^ 


'-\q\ 


271 


{qTp±{Q). 


M rn)\og{l/e) 

Cp,s = . - -r-, - P + 0 

47r \p\ 


(6.203) 


(6.204) 


There exists sfi dependence in the result (6.204) via n_s(0) = 1 — ns(0) 1/2 — 

(l/4)s/3/i, which can be easily understood as follows. The soft-fermion contribution to 
the damping rate comes from the process where a hard fermion (of type s) becomes 
a soft-fermion (of the same type s) by emitting a hard photon with almost the same 
momentum. The rate is proportional to {ubUpI) + 1)(1 — ^ 5 ( 0 )) where (1 — 795 ( 0 )) is 
the Pauli blocking factor of the final soft-fermion state, where we can put zero for soft- 
momentum at leading order in coupling. A similar process is where a hard fermion (of type 
s) meets with a soft-antifermion (of the type —s) to become a hard photon, and this rate is 
proportional to ( 79 b(|p|)-|- 1 ) 79 _ 5 ( 0 ) where 79 _ 5 ( 0 ) is the number density of soft-antifermion. 
The time reversed processes also add up to the damping rate, which is a property of 
fermionic case. These are each proportional to 79 b(|p|) 79 s( 0 ) and 79 b(|p|)(1 — 79 - 5 ( 6 )). 
Using 79 - 5 ( 0 ) - 1 - 795 ( 0 ) = 1, the total sum can be found to be 

(79 b(|p|) -M )(1 - 795(0)) (79 b(|p|) -h 1)79-5(0) -h 79^(1^1)795(0) 79 b(|p |)(1 - 79-5(0)) 

= 2 (r9B(|p|)+79-5(0)), (6.205) 


which nicely explains our result (6.204). 

II We point out that what is called rrij in Ref. [53] is in fact plasmino frequency which is equal to 
m^/2 in terms of asymptotic thermal mass m^. 
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Appendix 3: Expression of Fg^p, q; k) 

The function -Fs(p, qr; k) is given by 

q] k) = — _ __ _ _ _, (6.206) 

4 ll + p ■ p + kj \p\\p + q\\p + k\\p + q + k\ 

where 

^ = {\P\\P + q\ + \p\'^ + p ■ q) {\p + k\\p + k + q\ + \p + k\‘^ + {p + k) ■ q) 

+ {\P\\P + k\ + IpI^ +p - k) {\p + q\\p + q + k\ + \p + gp + (p + q) ■ k) 

- {\p\\P + q + k\ - p- {p + q + k)) {\p + q\\p + k\ - {p + q) ■ {p + k)) 

+ is {\p\ + \p + q\ + \p + k\ + \p + q + k\) e^^’'p'^q^k’’. (6.207) 

For the reasons mentioned in the main text, we are only interested in this quantity to 
linear order in the external momentum k and to second order in the loop momentum q. 
To this order the function Fs{p, q] k) is given by 


F,{p, q-k)r^l + ^ {{p • qf - \q\^) 

+ {p ■ ^ - 2(p ■ qf) + (p ■ k){p ■ q) - Ssie^^^Piqjkiip ■ q)) . 

(6.208) 


Appendix 4: jj? correction to electric conductivity 


Our analysis in this work contains all the necessary ingredients to compute the full fj? 
correction to the usual P-even electric conductivity at leading log order. The electric 
conductivity is given from Xs(|p|) by 


e 


2 


a = — 


3 



Xs(|p|), 

pO=s|p| 


(6.209) 


where Xs(|p|) satishes the second order differential equation written in (5.89), 


Cp,sXs(|p|) = S- 


e^m|) log(l/e) 
4t[ 


m 


rXs(|p|) 


m 


+ n.(lp|)-))x'.(lpl)-^x;(|p|) 

( 6 . 210 ) 
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where the soft-fermion contribution to the damping rate is given by (6.204), 


, log(l/e) 

Qs = jz ^ — (^b(IpI) + ^-^(0)) 


47r 


\P\ 


( 6 . 211 ) 


To correctly take into account corrections, we need to restore full expressions for m 
and mj including corrections. 
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TT^ 


^ 

TT^ 


( 6 . 212 ) 


Also, we have to expand ns(|p|) in (6.210) and n_s(0) in (6.211) up to second order in p. 


«.(|p|) = no{\p\ - s/i) no(|p|) - spn'odpl) + -n^UpDfi^ , 


(6.213) 


and n_s(0) ^ 1/2 — sl^jd/A + 0{n^), where no{x) = l/(e^* -1- 1). The resulting Xs(|p|) 
should be found up to order as 

Xs{\p\) = sX(o)(|p|) +hX(i)(|p|) + Vx( 2 )(|p|) + 0(p^), (6.214) 

where (x(o), X(i)) X( 2 )) can be obtained from (6.210) by solving it order by order in fi. 
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